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In this paper, we introduce nonlinear large- 
signal scattering (^) parameters, a new 
type of frequency-domain mapping that 
relates incident and reflected signals. We 
present a general form of nonlinear large- 
signal ^-parameters and show that they 
reduce to classic iS-parameters in the 
absence of nonlinearities. Nonlinear large- 
signal impedance (^ and admittance (^ 
parameters are also introduced, and equa- 
tions relating the different representations 
are derived. We illustrate how nonlinear 
large-signal ^-parameters can be used as a 
tool in the design process of a nonlinear 
circuit, specifically a single-diode 1 GHz 
frequency-doubler. For the case where a 
nonlinear model is not readily available, 
we developed a method of extracting non- 
linear large-signal ^-parameters obtained 
with artificial neural network models 
trained with multiple measurements made 



by a nonlinear vector network analyzer 
equipped with two sources. Finally, non- 
linear large-signal ^-parameters are com- 
pared to another form of nonlinear map- 
ping, known as nonlinear scattering func- 
tions. The nonlinear large-signal ^param- 
eters are shown to be more general. 
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1. Introduction 

Vector network analyzers (VNAs) are one of the 
most versatile instruments available for RF and 
microwave measurements. They are used to measure 
complex scattering parameters (^S-parameters) of linear 
devices or circuits. RF engineers use them to verify 
their designs, confirm proper performance, and diag- 
nose failures. A VNA works by exciting a linear device 
under test (DUT) with a series of sine wave signals, one 



frequency at a time, and detecting the response of the 
DUT at its signal ports. Since the DUT is linear, the 
input and output signal frequencies are the same as the 
source; these signals can be described by complex 
numbers that account for the signals' amplitudes and 
phases. The input-output relationships are described by 
ratios of complex numbers, known as ^S-parameters. For 
a two-port network, four ^S-parameters completely 
describe the behavior of a linear DUT when excited by 
a sine wave at a particular frequency. Although the 
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measurement of ^S-parameters by VNAs is invaluable to 
the microwave designer for modeling and measuring 
linear circuits, these measurements are oftentimes inad- 
equate for nonlinear circuits operating at large-signal 
conditions, since nonlinearities transfer energy from 
the stimulus frequency to products at new frequencies. 
Thus, conventional linear network analysis, which 
relies on the assumption of superposition, must be 
replaced by a more general type of analysis, which we 
refer to as nonlinear network analysis. 

Nonlinear network analysis involves characterizing a 
nonlinear device under realistic, large-signal operating 
conditions. To do this, complex traveling waves (rather 
than ratios) are measured at the ports of a DUT not only 
at the stimulus frequency (or frequencies), but also at 
other frequencies where energy may be created. 
Assuming the input signals are sine-waves and the 
DUT exhibits neither sub-harmonic nor chaotic behav- 
ior, the input and output signals will be combinations of 
sine-wave signals, caused by the nonlinearity of the 
DUT in conjunction with impedance mismatches 
between the measuring system and the DUT. If a single 
excitation frequency is present, new frequency compo- 
nents will appear at harmonics of the excitation fre- 
quency, and if multiple excitation frequencies are pres- 
ent, new frequency components will appear at the inter- 
modulation products as well as at harmonics of each of 
the excitation frequencies. In practice, there will be a 
limited number of significant harmonics and intermod- 
ulation products. The set of frequencies at which ener- 
gy is present and must be measured is known as the fre- 
quency grid. 

A class of instruments known as nonlinear vector 
network analyzers (NVNA) are capable of providing 
accurate waveform vectors by acquiring and correcting 
the magnitude and phase relationships between the fun- 
damental and harmonic components in the periodic sig- 
nals [1-5]. An NVNA excites a nonlinear DUT with one 
or more sine wave signals and detects the response of 
the DUT at its signal ports. Assuming the DUT does not 
exhibit any sub-harmonic or chaotic behavior, the input 
and output signals will be combinations of sine wave 
signals due to the nonlinearity of the DUT in conjunc- 
tion with mismatches between the system and the DUT. 
With these facts in mind, the major difference between 
a linear VNA and an NVNA is that a VNA measures 
ratios between input and output waves one frequency at 
a time while an NVNA measures the actual input and 
output waves simultaneously over a broad band of fre- 
quencies. 

Even though ^'-parameters cannot adequately repre- 
sent nonlinear circuits, some type of parameters relat- 



ing incident and reflected signals are beneficial so that 
the designers can "see" application-specific engineer- 
ing figures of merit that are similar to what they are 
accustomed to. In first part of this paper, we propose 
definitions of such ratios that we refer to as nonlinear 
large-signal scattering (^) parameters. We also intro- 
duce nonlinear large-signal impedance (^ and admit- 
tance (<^ parameters, and present equations relating 
the different representations. Next, we make two sim- 
plifications when considering the cases of a one-port 
network with a single-tone excitation and a two-port 
network with a single-tone excitation. 

For existing nonlinear models, we can readily gener- 
ate nonlinear large-signal ^-parameters by performing 
a harmonic balance simulation. For devices, with no 
model available, we can extract these parameters from 
artificial neural network (ANN) models that are trained 
with multiple frequency-domain measurements made 
on a nonlinear DUT with an NVNA. To illustrate appli- 
cations and generation of nonlinear large-signal ^ 
parameters, we present two examples. First, we illus- 
trate how nonlinear large-signal ^parameters can be 
used as a tool in the process of designing a simple non- 
linear circuit, specifically a single-diode 1 GHz fre- 
quency-doubler circuit. And secondly, we describe a 
method for generating nonlinear large-signal ^param- 
eters based upon ANN models trained on frequency- 
domain data measured using an NVNA. We compare a 
diode circuit model, generated using this method, to a 
harmonic balance simulation of a commercial device 
model. 

Finally, we compare our nonlinear large-signal ^- 
parameters to another form of nonlinear mapping, 
known as nonlinear scattering functions [6-7]. 
Specifically, we show that the two formulations are not 
equivalent. Nonlinear large-signal ^-parameters are 
more general than the nonlinear scattering functions, 
which are useful in approximating a specific class of 
nonlinearity in a more compact form. 



2. Nonlinear Large-Signal Scattering 
Parameters 

In this section, we introduce the concept of nonlinear 
large-signal scattering parameters. Like commonly 
used linear ^S-parameters, nonlinear large-signal scatter- 
ing (^) parameters can also be expressed as ratios of 
incident and reflected wave variables. However, unlike 
linear ^S-parameters, nonlinear large-signal ^-parame- 
ters depend upon the signal magnitude and must 
account for the harmonic content of the input and out- 
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put signals since energy can be transferred to other fre- 
quencies in a nonlinear device. 

After presenting the general form of nonlinear large- 
signal ^-parameters, we also introduce nonlinear 
large-signal impedance (M) and admittance (^ param- 
eters, and present equations for relating the different 
representations. Next, we make two simplifications in 
which we consider the cases of a one-port network with 
a single-tone excitation and a two-port network with a 
single-tone excitation. 

2.1 General Form 



For each nonlinear large-signal scattering parameter 
S>ijki the index / refers to the port number of the b wave, 
the indexy refers to the port number of the a wave, k is 
the harmonic index of the b wave, and / is the harmon- 
ic index of the a wave. The vectors a^ and b. are 
(M=3)-element vectors given by 



(5) 
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Equation (3) can be expanded as follows 



Consider an A^-port network. Normalized wave vari- 
ables aji and bji at theyth port and /th harmonic are pro- 
portional to the incoming and outgoing waves, respec- 
tively, and may be defined in terms of the voltages 
associated with these waves as follows: 



n 



■,b„ 



v: 



(1) 



where V^ji and V~ji represent voltages associated with 
the incoming and outgoing waves in the transmission 
lines connected to theyth port and containing frequen- 
cies of the /th harmonic; Z^j represents the characteris- 
tic impedance of the line at theyth port. 

The nonlinear large-signal scattering matrix ^ of the 
network expresses the relationship between a's and Z?'s 
at various ports and harmonics through the matrix 
equation 
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(6) 



Note that in each of the four sub-matrices, the diagonal 
elements contain the same-frequency scattering param- 
eters, the upper right elements contain the frequency 
down-conversion scattering parameters, and the lower 
left elements contain the frequency up-conversion scat- 
tering parameters. If the device under consideration 
contains no nonlinearities (i.e., no power is transferred 
to other frequencies), then Eq. (6) reduces to 



b=Ba, 



(2) 



where b and a are (A'^ x A/)-element column vectors. 
Here A'^ refers to the number of ports and M refers to the 
number of harmonics being considered. Matrix ^ is an 
(A'^x A/)^-element square matrix. We assume all a's and 
Z?'s are phase referenced to an to enforce time invari- 
ance [8]. 

As an example, consider a two-port network with 3 
harmonics; Eq. (2) then becomes 
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(7) 



which is the matrix representation for the well-known 
linear yS'-parameters involving three excitation frequen- 
cies. 

2.2 Nonlinear Large-Signal Impedance 
Parameters 

Rather than expressing the relationship between a's 
and Z?'s in terms of a nonlinear large-signal scattering 
matrix ^, we can alternatively express the relationship 
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between voltages (Fs) and currents (/'s) in terms of a 
nonlinear large-signal impedance matrix ^ as follows 



K=^/, 



(8) 



where V and / are (A^xM)-element column vectors. 
Once again A^ refers to the number of ports and M refers 
to the number of harmonics being considered. Z is an 
(A^xA/)^-element square matrix. 

For a two-port network with 3 harmonics, Eq. (8) 
becomes 






(9) 



where 



[^,j= 



-yll 



"y2i 



-ij3l 



"yi2 



^ij 22 



"y32 



^yi3 



^ij 23 



-ij33 



(10) 



^*=^;+^;=V^(%+^.x 



(13) 



where the subscripts refer to the /th port and the Ath har- 
monic. We can similarly express Iji as 

ij<=i^+ii~(!^-^)=^(^<-M (14) 



where the subscripts refer to the y'th port and at the /th 
harmonic. 

For simplicity, we will assume for now that the net- 
work under consideration consists of two ports. Later, 
we can easily generalize the equations relating the ^ 
and Z matrices for any A^-port network. If we allow the 
two transmission lines or waveguides connecting the 
two ports to have different characteristic impedances, 
Z^i and 7^2? Eq. (14) can be expressed in matrix form as 
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For each nonlinear large-signal impedance parameter 
Z^^i, the index / refers to the port number of the voltage 
F, the indexy refers to the port number of the current /, 
k is the harmonic index of V, and / is the harmonic 
index of /. The vectors V. and I. are (M=3)-element 
vectors given by 
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Equation (9) can be expanded to 
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2.3 Relating ^ and M Matrices 



The ^ and Z matrices can be expressed in terms of 
one another, if we know how a and b relate to V and /. 
From Eq. (1), we can express F^^ in terms of a^; and bn, 
as follows: 



where [U] is the identity matrix. Equation (9) can be 
expressed as 
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is the normalized impedance matrix. Equation (18) can 
be rewritten as 
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and Eq.(3) can be rewritten as 
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Combining Eqs. (20) and (21) allows us to solve for i 
in terms of ^. 
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If Z^i = Z„2? Eq. (22) reduces to 
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Alternatively, we can combine Eqs. (20) and (21) to 
solve for M in terms of d: 
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If Z^i = 2^2? Eq. (24) reduces to 



(24) 



[^11] [^2]" 

L iij L-12J ^25) 

[^21] [^22]. 



2.4 Nonlinear Large-Signal Admittance 
Parameters 

We can also express the relationship between volt- 
ages (Vs) and currents (Fs) in terms of a nonlinear 
large-signal admittance matrix i^, as follows 
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where ,^ is an (iVxA/)^-element square matrix. For a 
two-port network with three harmonics, for example, 
Eq. (26) becomes 
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For each nonlinear large-signal admittance parameter 
^yj,i, the index / refers to the port number of the current 
/, the indexy refers to the port number of the voltage F, 
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k is the harmonic index of/, and / is the harmonic index 
of V. The vectors Vj and /. are, once again, (M=3)-ele- 
ment vectors, defined in Eq. (11). Equation (27) can be 
expanded as follows 
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2.5 Relating ^ and ^ Matrices 



The ^ and ff matrices can also be expressed in 
terms of one another, using Eqs. (13) and (14) which 
show how a and b relate to V and /. 

For simplicity, we will again assume the network 
under consideration consists of two ports. If we allow 
the two transmission lines or waveguides connecting 
the two ports to have different characteristic imped- 
ances Z^i and Z„2? Eq. (14) can be expressed in matrix 
form as 
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where [U] is the identity matrix. Equation (27) can be 
expressed as 
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Combining Eqs. (30) and (31) gives 
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is the normalized admittance matrix. Equation (33) can 
be rewritten as 
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and Eq.(3) can be rewritten as 
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Combining Eqs. (35) and (36) allows us to solve for ^ 
in terms of ,^: 
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If Zo, = Zo2, Eq. (37) reduces to: 
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Alternatively, we can combine Eqs. (35) and (36) to 
solve for ^ in terms of ^: 
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IfZ^y = 2^2? Eq. (39) reduces to 
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2.6 One-Port Network With Single-Tone 
Excitation 

For a one-port network with a single-tone excitation 
at the fundamental frequency, we can extract a reflec- 
tion coefficient given by 
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The limitation imposed on the equation is that all other 
incident waves other than a^ equal zero. Instead of sim- 
ply taking the ratio ofbu, to a^, we reference the phase 
ofbik to that of an. To do this, we must subtract k times 
the phase of an from b^k [8]. 

For a one-port network with a single-tone excitation 
at the fundamental frequency, we can show that the 
equation relating ^ and Z reduces to the same well- 
known equation for the linear case if we assume that no 
energy is redistributed into the form of frequency 
down-conversion. To illustrate this, we will consider 
only M=3 harmonics, for the sake of simplicity. 
Equation (6) reduces to 
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for a one-port network with a single-tone excitation a^. 
This matrix can be rewritten as a set of three equations: 

^11=^1111^1 ; ki =^121^1 ; ^3 =^131^1- (43) 

Likewise, Eq. (12) reduces to 
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where the voltage V^ at the first harmonic can be 
expressed as 
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Combining Eqs. (45) and (46) gives 

V^(«ii+^ii)= 

[^llll(^l-Al)-^112^2 -^113^3 ]■ (47) 



Substituting Eq. (43) into Eq. (47) and solving for 
^1111 gives 

^ _ -^ol + ^1111 )+ ^1112 ^121 +^113^131 ^j^r.^ 

^1111 - /I O, \ ' ^ ^ 

(l-^llll) 

If no energy is redistributed into the form of frequency 
down-conversion (i.e., ^1112 = ^1113 = 0), then Eq. (48) 
reduces to the same equation as in the linear case: 
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7=7 
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A similar derivation can be performed to show that 
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Once again, if no energy is transferred to frequency 
down-conversion (i.e., #in2 = #1113 = 0), then Eq. (50) 
reduces to the same equation as in the hnear case: 



^^1 = 



1 _ 1 O-^n) 
^11 ^.1 + ^11)' 



(51) 



2.7 



Two-Port Network With Single-Tone 
Excitation 



For a two-port network excited at port 1 by a single- 
tone excitation at the fundamental frequency, we can 
extract an input reflection coefficient given by 
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As with Eq. (41), instead of simply taking the ratio of 
bik to an, we phase reference to a^. To do this we must 
subtract ^ times the phase of an from bn,. The limitation 
once again imposed on the equation is that all other 
incident waves other than a^ equal zero. 

Another valuable parameter, the forward transmis- 
sion coefficient, is similarly extracted as follows 



KKt^^-ktJ 



DUT, one of the equations in the matrix defined by Eq. 
(6) reduces to 



^2k -^2lkAl '^^llkk^k 



(54) 



If we solve Eq. (54) for ^22^5 we obtain 



1^1 



(55) 

^2k ^2k 

In Eq. (55), the output reflection coefficient ^22kk obvi- 
ously cannot be determined by simply taking the ratio 
ofb2k to a2k, since the ratio also depends on a^ through 
^2iM- When a2k is small, we can generate another signal 
Aa2jt that is offset slightly from the frequency of interest 
ff, by A/^. Eq. (54) then becomes 

hk -^^b^k = ^mi^i +^22^ {^k +^% X (^^) 

where Aa2;s, « a2;s, and ^22^: remains constant over this 
frequency range. Subtracting Eq. (54) from Eq. (56) 
gives 



AZ>.. 



k^^2k^ 



(57) 



which does not depend on ^2ui- If we solve Eq. (57) 
for ^22kk, we obtain 



_A^ 



Aa,, 



Large a^, Small Aa2}, 



(58) 



Equation (58) is a quasi-linear approximation of the 
output reflection coefficient under normal operating 
conditions, and is consistent with the definition of "Hot 
^'22," which has been used to measure the degree of mis- 
match at the output port of a power amplifier at its exci- 
tation frequency. 



a^^^ =0 for all m, n ^^\)A(n^l)\ 2.8 Summary of Sec. 2 



(53) 

This parameter provides a value of the gain or loss 
through a device either at the fundamental frequency or 
converted to a higher harmonic frequency. 

In addition to the previous two parameters, given in 
Eqs. (52) and (53), an output reflection coefficient can 
also be useful when trying to determine the output 
matching network. If a nonlinear DUT is operating 
under its normal drive condition (an at some constant 
signal level), and a second source, excited by a small- 
signal tone at frequency f^, is placed at port 2 of the 



In this section, we presented the general form of non- 
linear large-signal ^-parameters. Unlike linear 5*- 
parameters, nonlinear large-signal ^-parameters 
depend upon the signal magnitude and must take into 
account the harmonic content of the input and output 
signals, since energy can be transferred to other fre- 
quencies in a nonlinear device. We also introduced non- 
linear large-signal impedance (^ and admittance (#) 
parameters, and presented equations for relating the 
different representations. Next, we made two simplifi- 
cations, considering the cases of a one-port network 
with a single-tone excitation and a two-port network 
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with a single-tone excitation. For the one-port case with 
a single-tone excitation at the fundamental frequency, 
we showed that the equation relating ^ and M reduces 
to the same well-known equation for the linear case if 
we assume that no energy is transferred to frequency 
down-conversion. For the two-port case excited at port 
1 by a single-tone excitation at the fundamental fre- 
quency we extracted an input reflection coefficient 
^iiM? a forward transmission coefficient ^2ui? and a 
quasi-linear output reflection coefficient ^2m- 



3, Using Nonlinear Large-Signal ^ 
Parameters to Design a Diode 
Frequency-Doubler Circuit With a 
Harmonic-Balance Simulator 

Resistive frequency doublers operate on the princi- 
ple that a sinusoidal waveform is distorted by the non- 
linear //F characteristic of a Schottky-barrier diode [9]. 
This distortion causes power to be generated at higher- 
harmonic frequencies. The design of such doublers 
involves separating the input and output signals by fil- 
ters and determining the optimum input and output 
matching circuits, as illustrated in Fig. 1. 



Z^ Input 
~H Matcliing 
('^\ Network 





-^^4- 




Irput 
Filtering 

Network 




OuLpul 
Filtering 
Network 



Output 
Matcliing 
Network 



Zr 



Fig. 1. Block diagram of a single-diode resistive doubler. 

Although single-diode resistive doublers are not very 
efficient (analysis predicts a conversion loss of at least 
9 dB [10]), we chose this circuit because it is simple 
enough to clearly illustrate how nonlinear large-signal 
^-parameters can be used as a design tool. 

In the following sections, we describe the various 
steps involved in designing a single-diode 1 GHz fre- 
quency-doubler circuit. Since we are using a simulator, 
we can force the stimulus to consist of only |aii|, with 
all other a^^ terms equal to zero, where m and n are pos- 
itive integers such that m ^ 1 and n^ \. (In practice, 
this condition can never be completely realized in a 
measurement environment.) With only an a^ compo- 
nent present, we need only consider the parameters 
^lui (Eq. 52), which is a measure of the large-signal 
input match at the Ath harmonic, as well as the param- 



eter ^2i;ti (Eq. 53), a measure of the large-signal con- 
version loss or gain at the Ath harmonic, plus the quasi- 
linear ^2222 (Eq. 58) to determine the output matching 
network at the second harmonic. Figure 2 illustrates the 
setups required for determining these parameters. 
Determining ^2222 requires a second source at port 2 at 
a frequency slightly offset from 6;^. 



7 


Nonlinear 
DOT 










Zl 









00 




Fig 2. Nonlinear large-signal ^parameters used to characterize a 
two-port device excited by a single-tone signal at port 1. 

In the first step, we perform a simulation on the 
diode alone and use ^2121 to determine the optimum 
bias condition for converting power from the funda- 
mental frequency to the second harmonic. Second, we 
add filtering networks to separate the input and output 
signals, and verify their proper performance by looking 
at ^2111 and ^1121- Third, we make use of ^im to deter- 
mine the input matching network. Fourth, with the 
input matching network in place, we place a second 
source at port 2 and find the quasi-linear value of ^22225 
which allows us to determine the output matching net- 
work. Fifth, we use the optimization feature of the sim- 
ulator to minimize ^im by varying the line lengths of 
the input and output matching circuits. And finally, 
sixth, we add 4 GHz and 6 GHz filters at the output 
(and re-determine the proper input and output matching 
circuits) in order to reduce the values of ^2141 and ^2161? 
which in turn increases the value of ^2121 and cleans up 
the output waveform. 

3.1 Diode Only 

In this example, we use a compact model to simulate 
a commercial Schottky-barrier diode. The model 
includes a series resistance T^g of 14 Q, a junction 
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capacitance at zero voltage Cjq of 0.08 pF, and a reverse 3.3 

saturation current /g of 3 x 10"^^ A. 

First, we perform a harmonic-balance simulation on 
the diode, sweeping the bias voltage to determine 
which condition gives the highest value of ^2121 for 
an = 1.0 V. Note that in all simulations we set the gen- 
erator impedance Zq and the load impedance Zl to 
50 ^. After sweeping the voltage, we determine that the 
optimum forward bias is +0.48 V. 

3.2 Diode With 1 GHz and 2 GHz Filters 

With a stimulus of an = 1.0 V and a forward bias of 
+0.48 V, we add filtering networks to separate the input 
and output signals. On the input side, we place a 2 GHz, 
A/4 (A/8 at 1 GHz) open-circuited stub. This creates an 
RF short at 2 GHz, preventing the output power gener- 
ated in the diode from traveling backward. On the out- 
put side, we place a 1 GHz, A/4 open-circuited stub. 
This creates an RF short at 1 GHz, preventing any sig- 
nal at 1 GHz from traveling forward. 

Table 1 lists the simulated values for ^nn - ^u^i, 
^2111 - ^2161? ^2 and <§2^<§ for each of the design 
stages, where # is the expanded power gain and #2 is 
the expanded power gain confined to the second har- 
monic, as defined in [11]. With the 1 GHz and 2 GHz 
filters in place, we see that the value of l^i^il decreas- 
es from 0.170 to 1.3 x 10"^ the value of |^2iiil decreas- 
es from 0.536 to 3.3 x 10"^ and (§2 increases from 
-14.16 dB to -9.73 dB. 



Diode With 1 GHz and 2 GHz Filters and 
Input Matching 



Once the filters are placed in the circuit, we make use 
of the complex-valued ^im to design the input match- 
ing network with the well-known single, open-circuited 
stub technique. This is possible, assuming that no ener- 
gy is transferred to frequency down-conversion, as dis- 
cussed in Sec. 2.6. We see in Table 1 that l^iml reduces 
from 0.569 without the input matching network to 
9.4 X 10"^ with the input matching network in place. 
Likewise, (§2 increases from -9.73 dB to -9.69 dB. 

3.4 Diode With 1 GHz and 2 GHz Filters, Plus 
Input and Output Matching 

Whereas our input matching network is designed for 
1 GHz, our output matching network must be designed 
for 2 GHz. While the circuit is operating under its nor- 
mal drive condition (an = 1.0 V and a forward bias of 
+0.48 V) we place a second source at port 2, excited by 
a small-signal tone (Aa22 = 0.01 V) at a frequency off- 
set of 10 kHz from the desired 2 GHz, to give us the 
quasi-linear value of ^22225 which allows us to deter- 
mine the output matching network. We make use of 
^2222 to design the output matching network with the 
well-known single, open-circuited stub technique. We 
see in Table 1 that with the output matching network in 
place, the value of 1^2121 1 is only marginally increased 
from 0.326 to 0.328. This is because the value of ^2222 



Table 1. Simulated values for ^ 
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^2> and ^2^^ for each of the design stages of the diode frequency doubler 
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is relatively low, which means the output is already 
almost matched to 50 Q. We also note that <^2 increas- 
es from -9.69 dB to -9.65 dB. 

3.5 Diode With 1 GHz and 2 GHz Filters, Plus 
Optimized Input and Output Matching 

With the filters and matching networks in place, we 
use the optimization feature of the simulator to mini- 
mize ^1111 by varying the lengths of the lines in the 
input and output matching circuits. Doing this decreas- 
es the value of |^uu| from 8.7 x 10"^ to 6.0 x 10"^ while 
increasing the value of \^2m\ from 0.328 to 0.331 and 
(§2 from -9.65 dB to -9.60 dB. 
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Filters 
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Fig. 3. Final design of the single-diode resistive frequency doubler. 
Electrical lengths shown are all at 1 GHz. 

3.7 Summary of Sec. 3 



3.6 Diode With (1, 2, 4, and 6) GHz Filters, Plus 
Optimized Input and Output Matching 

From Table 1, we see that at the output port, |^2iiil? 
I^2i3il5 and |^2i5il all have values less than or equal to 
4.0 X 10"^, but |^2i4il and |^2i6il have noticeably higher 
values (at least 2.9 x 10"^). 

In order to clean up the output waveform, we add 4 
GHz and 6 GHz filters, in the form of A/4 open-circuit- 
ed stubs, at the output. With these filters placed in the 
circuit, we re-determine the proper input and output 
matching conditions. After optimizing the circuit once 
again, the value of |^2i4il decreases from 4.0 x 10"^ to 
1.4x10"^ and the value of |^2i6il decreases from 
2.9 X 10"^ to 2.7 X 10"^ The addition of these filters, in 
turn, slightly increases |^2i2il froin 0.331 to 0.332 and 
<§2 from -9.60 dB to -9.56 dB. At this final design 
stage, the overall power gain is nearly -9.56 dB since 
the ratio <§2^<§= 0.999. The semi-empirical analysis of 
[10] predicts a maximum gain of -9 dB. Figure 3 illus- 
trates the final design of the single-diode resistive dou- 
bler circuit. And Fig. 4 shows the time-domain plots of 
ai and b2 for the final design of the simulated 1 GHz 
frequency-doubler circuit. 




We illustrated how nonlinear large-signal ^param- 
eters can be used as a tool in the design process of a sin- 
gle-diode 1 GHz frequency-doubler. Specifically, we 
used ^1111 to determine the input matching network, 
^2222 to determine the output matching network, and 
^iijti, ^2iH (for ^ = 1 to 6), and <§2 to quantify the per- 
formance of the circuit at each stage. 

By the final stage of the design, we had created a 
doubler with an overall power gain of -9. 56 dB, not far 
from the maximum possible predicted value of -9 dB. 



4. Determining Nonlinear Large-Signal 
^Parameters from Artificial Neural 
Network Models Trained With 
Measurement Data 

Although nonlinear large-signal ^-parameters can 
be easily determined for an existing model in a com- 
mercial harmonic balance simulator by forcing all a's 
other than an to zero, they cannot be determined direct- 
ly from measurements. With currently available 
NVNAs, the nonlinear DUT, in conjunction with the 
impedance mismatches and harmonics from the system 
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Fig. 4. Time-domain plots ofa^ and ^2 for tlie simulated 1 GHz frequency-doubler circuit. 
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make it impossible to set all a's other than a^ (assum- 
ing port 1 excitation) to zero. In order to overcome this 
obstacle, we propose a method [12] that makes use of 
multiple measurements of a DUT using a second source 
with isolators, as shown in Fig. 5. This measurement 
set-up is similar to that introduced by Verspecht et al. 
[6-7] to generate "nonlinear scattering functions." As a 
side note, we compare and contrast the "nonlinear scat- 
tering functions" with our definitions of nonlinear 
large-signal scattering parameters in the Appendix. 



Precidon ADC 



Droadtand DoTSTiconverter 



DiasT 



Kk ^ik 



DLTT 




o^ 






Fig. 5. Block diagram of a nonlinear vector network analyzer 
equipped with a second source and isolators. 

4.1 Methodology 

To illustrate our technique of generating nonlinear 
large-signal ^-parameters, let us consider the case 
where a DUT is excited at port 1 by a single-tone sig- 
nal at frequency /i and signal level |aii|. Utilizing a sec- 
ond source, we take multiple measurements of a non- 
linear circuit for different values of a^„ [{mi^\)A{n-i^\)'\. 
We then use these data to develop an artificial neural 
network (ANN) model that maps values of a's to Z?'s, as 
shown in Fig. 6. Once the ANN model is trained and 
verified, the nonlinear large-signal ^-parameters are 
obtained by interpolating Z?'s from the measured results 
for nonzero values of a^„ [{mi^\)A{ni^\)'\ to the desired 
values for a^^ [{mi^\)A{ni^\)'\ equal to zero, as shown in 
Fig. 7. Alternatively other conditions maybe called for, 
where a^„^0 depending on the desired application-spe- 
cific figure of merit. 

One popular type of ANN architecture, which is used 
in our work, is a feed-forward, three-layer perceptron 
structure (MLP3) consisting of an input layer, a hidden 
layer, and an output layer [13]. The hidden layer allows 
for complex models of input-output relationships. 
ANNs learn relationships among sets of input-output 
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Fig. 6. An ANN model that maps real and imaginary values of a's to 
^'s for different real and imaginary values of a^^ [{m^V)A.{n^\)\ 
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Fig. 7. An ANN model that interpolates ^'s from the measured 
results for nonzero values of a^^ [(m^l)A(n^iy] to the desired values 
for fl^„ [(m?^l)A(?/?^l)] equal to zero. Outputs of the ANN model 
yield values of ^^ij^i. 



data that are characteristic of the device or system 
under consideration. After the input vectors are present- 
ed to the input neurons and output vectors are comput- 
ed, the ANN outputs are compared to the desired out- 
puts and errors are calculated. Error derivatives are 
then calculated and summed for each weight until all of 
the training sets have been presented to the network. 
The error derivatives are used to update the weights for 
the neurons, and training continues until the errors 
become no greater than prescribed values. In our study, 
we have utilized software developed by Zhang et al. 
[14] to construct our ANN models. 

To test our method of generating nonlinear large-sig- 
nal ^-parameters, we fabricated a wafer-level test cir- 
cuit using a Schottky diode in a series configuration, as 
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shown in Fig. 8. The two-port diode circuit was fabri- 
cated on an alumina substrate by bonding a beam-lead 
diode package to the gold metalization layer with silver 
epoxy The diode was located in the middle of the 
coplanar waveguide (CPW) transmission lines, with 
short lines connecting the diode to probe pads at both 
ports. We measured the test circuit on an NVNA using 
an on-wafer VNA line-reflect-reflect-match (LRRM) 
calibration, along with signal amplitude and phase cal- 
ibrations. This process places the reference plane at the 
tips of the wafer probes used to connect with the CPW 
leads. 



J- 



h- 



Fig. 8. Schottky diode in a series configuration located in the mid- 
dle of a CPW transmission line. (White area is metal.) 

For all measurements, the first source, located at port 
1 , used a sine-wave excitation of frequency 900 MHz 
and magnitude |aii|==0.178 V (-5 dBm in a 50 ^ envi- 
ronment) at the probe tips. The second source was con- 
nected to port 2 and used a sine-wave excitation of fre- 
quency 900 MHz and |a2i|-0.178 V. The diode was for- 
ward-biased to +0.2 V through the probe tips. In order 
to obtain the nonlinear large-signal ^-parameters, ^^i 
and ^2iM? the excitation from source 1 was held con- 
stant, while the phase of source 2 was randomly 
changed for 500 different measurements that varied 
slightly in magnitude. Figure 9 plots the resulting 
measurements of a2i in the complex plane. The nonlin- 
earities in the test circuit, along with impedance mis- 
matches, created other input components at higher har- 
monics, as shown in Figs. 10-13 for the second and 
third harmonics (a^, ciu, cin^ and a^^). These variations 
in a^j allowed us to create an ANN model that could be 
used to interpolate Z?'s from the measured results for 
nonzero values of a^„ [(m^l)A(n^l)], as shown in Figs. 
14 and 15 for b^y and Z?2i, to the desired values for 
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Fig. 9. Five hundred measurements of 021 iii the complex plane with 
the excitation from source 1 held constant and the output from source 
2 set to random phases with constant amplitude. 
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Fig. 10. Five hundred measurements of 0^2 in the complex plane 
with the excitation from source 1 held constant and the output from 
source 2 set to random phases with constant amplitude. 

<^mn [(m^l)A(n^l)] equal to zero, or alternatively anoth- 
er desired device condition. 

4.2 Sensitivity Analysis of ANN Models 

Data from the 500 measurements were used to devel- 
op two ANN models, one for mapping values from the 
first five harmonics of ai and a2 (an, ai2, . . ., a^^, a2i, a22, 
. . ., a25) to the first five harmonics of Z?i (Z?ii, Z?i2, . . ., by^), 
and the other for mapping values from the first five har- 
monics of «! and a2 to the first five harmonics of Z?2 (^21? 
^225 ---. ^25)- We performed a sensitivity analysis to 
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Fig. 11. Five hundred measurements of a^^ in the complex plane 
with the excitation from source 1 held constant and the output from 
source 2 set to random phases with constant amplitude. 
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Fig 12. Five hundred measurements of ^22 in the complex plane with 
the excitation from source 1 held constant and the output from source 
2 set to random phases with constant amplitude. 



Fig. 13. Five hundred measurements of ^23 in the complex plane 
with the excitation from source 1 held constant and the output from 
source 2 set to random phases with constant amplitude. 
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Fig 14. Five hundred measurements of /^n in the complex plane with 
the excitation from source 1 held constant and the output from source 
2 set to random phases with constant amplitude. 



determine how many training points, testing points, and 
hidden neurons are required to adequately train the two 
ANN models. Tables 2-4 summarize the results for the 
first model, where we map values from the first five 
harmonics of ai and aj to the first five harmonics of Z?i, 
and Tables 5-7 summarize the results for the second 
model, where we map values from the first five har- 
monics of ai and aj to the first five harmonics of Z?2. 

First, we varied the number of hidden neurons from 
1 to 20. All other parameters were held constant. 
Specifically, the 500 measurements points were divid- 
ed into 250 training points and 250 testing points, and 



we used the conjugate gradient method for training. 
Table 2 lists the average testing errors and correlation 
coefficients for the models that map ai and aj to Z?i, and 
Table 5 lists the average testing errors and correlation 
coefficients for the models that map ai and a2 to bj- 
Both mappings show similar trends. The average test- 
ing errors decreased with increasing numbers of hidden 
neurons until around 14 or 16, where the errors were 
minimized. For more than 1 6 hidden neurons, the trend 
reversed and the errors appeared to start increasing 
again. Figure 16 plots the average testing errors as a 
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Fig. 15. Five hundred measurements of ^21 ^^ the complex plane 
with the excitation from source 1 held constant and the output from 
source 2 set to random phases with constant amplitude. 



Table 3. Average testing errors and correlation coefficients as func- 
tions of the number of training points for ANN models trained to map 
values from the first five harmonics ofai and a2 to the first five har- 
monics of ^1. All models were developed using 14 hidden neurons 
and verified using 250 testing points 



Training 


Average testing 


Correlation 


points 


error (%) 


Coefficient 
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20.10 


0.96764 


10 


9.01 


0.99556 
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3.64 


0.99891 


50 


1.91 


0.99979 


125 


0.95 


0.99995 


250 


0.72 


0.99997 



Table 4. Average testing errors and correlation coefficients as func- 
tions of the number of testing points for ANN models trained to map 
values fiom the first five harmonics ofai and a2 to the first five har- 
monics of b]^. All models were developed using 250 training points 
and 14 hidden neurons 
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Fig. 16. Average testing errors as fiinctions of the number of hidden 
neurons for ANN models trained to map a^ and a2 to b^ and a^ and 
a2 to ^2- The models were developed using 250 training points and 
verified using 250 testing points. 



Table 2. Average testing errors and correlation coefficients as fiinc- 



Testing 


Average testing 


Correlation 


points 


error (%) 


Coefficient 


5 


0.80 


0.99998 


10 


0.74 


0.99997 


25 


0.68 


0.99998 


50 


0.68 


0.99998 


125 


0.72 


0.99997 


250 


0.72 


0.99997 



Table 5. Average testing errors and correlation coefficients as func- 
tions of the number of hidden neurons for ANN models trained to 
map values from the first five harmonics ofa^ and a2 to the first five 
harmonics of ^2- ^^ models were developed using 250 training 
points and verified using 250 testing points 



tions of the number of hidden neurons for ANN models trained to 


Hidden 


Average testing 
error (%) 


Correlation 


map values from the first five harmonics ofa^ and 


^2 to the first five 


neurons 


Coefficient 


harmonics 


of ^1. 
verified 


All models were developed using 250 training 
using 250 testing points 








points and 
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17.88 
13.22 


0.74320 










0.91161 


Hidden 




Average testing 


Correlation 
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6.48 


0.96659 


neurons 




error (%) 


Coefficient 


6 
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2.04 
1.43 


0.99893 










0.99951 
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16.86 


0.94814 
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0.99985 
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0.98896 


12 


0.82 


0.99989 


4 




4.56 


0.99715 
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0.78 


0.99989 


6 




1.66 


0.99971 


16 


0.73 


0.99992 
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1.15 


0.99989 
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0.78 


0.99988 


10 




1.08 


0.99991 


20 


0.99 


0.99983 
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0.80 
0.72 


0.99996 
0.99997 
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20 




0.70 


0.99997 
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Table 6. Average testing errors and correlation coefficients as func- 
tions of the number of training points for ANN models trained to map 
values from the first five harmonics ofa^ and aj to the first five har- 
monics of ^2- ^11 models were developed using 14 hidden neurons 
and verified using 250 testing points 



Training 
points 

5 

10 

25 

50 
125 
250 



Average testing 
error (%) 

27.08 
12.99 

3.72 

1.75 

1.09 

0.78 



Correlation 
Coefficient 

0.50237 
0.91962 
0.99628 
0.99940 
0.99978 
0.99989 



Table 7. Average testing errors and correlation coefficients as func- 
tions of the number of testing points for ANN models trained to map 
values from the first five harmonics ofa^ and a2 to the first five har- 
monics of ^2- All models were developed using 250 training points 
and 14 hidden neurons 



Testing 
points 



Average testing 
error (%) 



Correlation 
Coefficient 



5 

10 

25 

50 

125 

250 



0.87 
0.84 
0.81 
0.80 
0.81 
0.78 



0.99995 
0.99993 
0.99988 
0.99989 
0.99988 
0.99989 



function of the number of hidden neurons for both map- 
pings. 

Next, we varied the number of training points from 5 
to 250. All other parameters were held constant. The 
number of hidden neurons was set to 14 since we found 
that to be an ideal number from the previous analysis, 
and 250 testing points were used for verification. Table 
3 lists the average testing errors and correlation coeffi- 
cients for the models that map a^ and a2 to by, and Table 
6 lists the average testing errors and correlation coeffi- 
cients for the models that map a^ and a2 to b2. Once 
again, both mappings showed similar trends. The aver- 
age testing errors decreased for an increasing number 
of training points. However, as more and more training 
points were added, diminishing returns on the testing 
errors were evident. Figure 17 plots the average testing 
errors as a function of the number of training points for 
both mappings. 

Finally, we varied the number of testing points from 
5 to 250. All other parameters were held constant. The 
number of hidden neurons was once again set to 14, and 
the same 250 training points were used for model 
development. Table 4 lists the average testing errors 
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Fig. 17. Average testing errors as functions of the number of train- 
ing points for ANN models trained to map a^ and a2 to by and a^ and 
a2 to ^2- The models were developed using 14 hidden neurons and 
verified using 250 testing points. 

and correlation coefficients for the models that map a^ 
and a2 to b^, and Table 7 lists the average testing errors 
and correlation coefficients for the models that map a^ 
and a2 to Z)2. Both mappings showed that the average 
testing errors varied little with the number of testing 
points. Figure 18 plots the average testing errors as a 
function of the number of testing points for both map- 
pings. 
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Fig. 18. Average testing errors as functions of the number of testing 
points for ANN models trained to map a^ and a2 to b^ and a^ and a2 
to ^2- The models were developed using 14 hidden neurons and 250 
training points. 



4.3 Results and Comparison for Sec. 4 

Based on the resuhs of our sensitivity analysis, we 
decided to use 250 training points and 250 testing 
points to train and verify the two ANN models. We 
chose to use 14 hidden neurons for mapping values 
from the first five harmonics of ai and a2 to the first five 
harmonics of bi and 16 hidden neurons for mapping 
values from the first five harmonics of ai and a2 to the 
first five harmonics ofb2- The testing error was 0.72 % 
for the bi model and 0.73 % and for the b2 model, with 
respective correlation coefficients of 0.99997 and 
0.99992. 
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After the ANN models were developed, the nonlin- 
ear large-signal ^-parameters, ^i^.^ and ^2ijti (^= 1)2, 
. . ., 5), were obtained by interpolating bi^ and b2k from 
measured results for nonzero values of ai2, ^13, -.-, ^15 
and a2i, ^22, . . ., ^25 to the desired values for ai2, ^13, - - -, 
ai5 and a2i, ^22, -.-, ^25 equal to zero. Figure 19 shows 
the interpolated value of Z?ii (= ^im • a^) when ai2, ^13, 
..., ai5 and a2i, ^22, -.-, ^25 were set equal to zero, and 
Fig. 20 shows the interpolated value of Z?2i (= ^2111 ' ^11) 
when ai2, ^13, - - -, ^15 and a2i, ^22, - - -, ^25 were set equal 
to zero. 
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Fig, 19, The 250 measurements of b^^ used for training (circles). 
Values of ^1111 • «n were determined from the measurement-based 
ANN model (square) and the harmonic balance simulation using a 
compact model (triangle). 
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Fig, 20, The 250 measurements of b2i used for training (circles). 
Values of ^2111 ' ^11 ^^^^ determined from the measurement-based 
ANN model (square) and the harmonic balance simulation using a 
compact model (triangle). 



We compared our results to a compact model provid- 
ed by the manufacturer and simulated in commercial 
harmonic-balance software to get an independent check 
on our methodology. Our comparison was accom- 
plished by providing the simulator with the identical 
biasing conditions on the diode and a stimulus of the 
same magnitude used in the measurements for a^ and 
setting all other a's to zero. Providing the simulated cir- 
cuit with ail of the same magnitude as the measurement 
should give the same values of bu, and Z?2jt as the inter- 
polated values of Z?ijt (= ^nki ' ^11) and b2k (= ^2ih ' ^11) 
determined by the ANN models when ai2, ^13, ..., a^^ 
and a2i, ^22, . . ., ^25 are set equal to zero. Figures 19 and 
20 show that the simulated values Z?ii and Z?2i agree with 
those determined from the measurement-based ANN 
models. 

Quantitatively, the differences between the ANN and 
equivalent-circuit models are shown in Table 8. 

4.4 Summary of Sec. 4 

We described a method of extracting nonlinear large- 
signal ^parameters, using an NVNA equipped with 
isolators and a second source. First, we showed how 
multiple measurements of a nonlinear circuit could be 
used to train artificial neural networks. Then, we 
extracted the desired ^-parameters by interpolating the 
ANN models for all a's equal to zero other than an. We 
checked our approach by comparing our results to a 
compact model simulated in commercial harmonic-bal- 
ance software, and showed that the two methods agree 
well. 

We also performed a sensitivity analysis on the ANN 
networks, and discovered the following: (1) The aver- 
age testing error decreases for an increasing number of 
training points. However, as more and more training 
points are added, diminishing returns on the testing 
errors are evident. (2) As the number of hidden neurons 
are increased, the average testing error decreases until 
around 14 hidden neurons at which point more hidden 
neurons have no benefit and can actually lead to 
increases in testing error. (3) The number of testing 
points does not drastically affect the testing error. In 
fact, no more than 25 testing points are needed for the 
models tested. 



5. Overall Summary 

In this paper, we introduced nonlinear large-signal 
scattering parameters representing a new type of fre- 
quency-domain mapping that relates incident and 
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Table 8. Diflferences between the measurement-based, ANN-modeled results and the compact model simulated in commercial harmonic-balance 
software 



Quantity 


Diflference 


Diflference 


Quantity 


Diflference 


Difference 




(%) 


(dBV) 




(%) 


(dBV) 


^1111 


3.38 


-A4.5 


^2111 


3.95 


^3.2 


^1121 


1.23 


-53.3 


^2121 


7.15 


-38.0 


^1131 


3.29 


^4.8 


^2131 


5.93 


-39.6 


^1141 


0.40 


-63.1 


^2141 


0.72 


-57.9 


^1151 


1.67 


-50.6 


^2151 


0.85 


-56.5 



reflected signals. Unlike classical ^S-parameters, nonlin- 
ear large-signal ^-parameters take harmonic content 
into account and depend on the signal magnitudes. 
First, we presented a general form of nonlinear large- 
signal ^-parameters and showed that they reduce to 
classic iS-parameters in the absence of nonlinearities. 
We also introduced nonlinear large-signal impedance 
(^ and admittance (^) parameters, and presented 
equations that relate the different representations. Next, 
we considered two simplified cases of a one-port net- 
work and a two-port network, each with a single-tone 
excitation. For the one-port network, we showed that 
the equation relating S> and ^reduces to the same well- 
known equation for the linear case, assuming no power 
is transferred in the form of frequency down-conver- 
sion. For the two-port case, we extracted input reflec- 
tion coefficients and forward transmission coefficients, 
which can be useful for designing circuits such as 
amplifiers and frequency multipliers. In addition, we 
derived a quasi-linear approximation of the output 
reflection coefficient under normal operating condi- 
tions. These three two-port parameters allow a design- 
er to "see" application-specific engineering figures of 
merit that are similar to what he or she is accustomed to 
in the linear world. 

Next, we illustrated how nonlinear large-signal ^ 
parameters can be used as a tool in the design process 
of a single-diode 1 GHz frequency-doubler. 
Specifically, we used ^nii to determine the input 
matching network, ^2222 to determine the output match- 
ing network, and ^nku ^im (for ^= 1 to 6), and (§2 to 
quantify the performance of the circuit at each stage. 
By the final stage of the design, we had created a dou- 
bler with an overall power gain of -9. 56 dB, a value not 
far from the maximum possible predicted value of -9 
dB. 

For the case where a nonlinear model is not readily 
available, we described a method of extracting nonlin- 
ear large-signal ^-parameters, using an NVNA 
equipped with isolators and a second source. First, we 
showed how multiple measurements of a nonlinear cir- 



cuit could be used to train artificial neural networks. 
Then, we extracted the desired ^parameters by inter- 
polating the ANN models for all a's equal to zero other 
than an. We checked our approach by comparing our 
results to a compact model simulated in commercial 
harmonic-balance software, and showed that the two 
methods agree well. We also performed a sensitivity 
analysis on the ANN networks, and discovered the fol- 
lowing: (1) The average testing error decreases for an 
increasing number of training points. However, as more 
and more training points are added, diminishing returns 
on the testing errors are evident. (2) As the number of 
hidden neurons are increased, the average testing error 
decreases until around 14 hidden neurons, at which 
point more hidden neurons have no benefit and can 
actually lead to increases in testing error. (3) The num- 
ber of testing points does not drastically affect the test- 
ing error. In fact, no more than 25 testing points are 
needed for the models tested. 



6, Appendix A, Comparing Nonlinear 
Large-Signal ^Parameters With 
Nonlinear Scattering Functions 

Here, we compare the nonlinear large-signal ^ 
parameters, introduced in this paper, to another form of 
nonlinear mapping, known as nonlinear scattering 
functions, introduced by Verspecht [6-7]. 

For a two-port nonlinear device, excited by a single- 
tone signal, and assuming all harmonic signals are rel- 
atively small compared to the fundamental signals, 
Verspecht defines nonlinear scattering functions as 



K = K^ 



I ^. 



/=1,2 
j=2,...,M 



kpij 



Re(^^) + 



X H^,j\m{qj), {S9) 



i=\,2 
j=2,...,M 



where a^j and bf^, represent the wave variables propor- 
tional to the incoming and outgoing waves, respective- 
ly, and M refers to the number of harmonics being taken 
into account. F^^, Gj^y, and Hf^^y are functions of the fun- 
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damental components Re(aii), Re(a2i), and Im(a2i). The 
imaginary component of a^ is omitted, with the 
assumption that the wave variables are phase refer- 
enced such that the phase of an is set to zero. F^, Gj^ip 
and H^ij are assumed complex constants for a given 
bias and fundamental drive condition. Note that these 
three terms do not depend upon the higher harmonic 
signal levels. With the Uy wave variables split into real 
and imaginary components, G^y and H^y serve to map 
aij circles centered at zero to bj^ ellipses with variable 
axes also centered at zero, as shown in Fig. 21. The Fj^ 
terms translate the ellipses about the complex plane. 
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Fig. 21. Gj^y and H^y serve to map ay circles centered at zero to bj^ 
ellipses witii variable axes also centered at zero, neglecting F^ for 
illustrative purposes. 



For illustrative purposes, let us consider bjj, taking 
into account the first three harmonics. Doing this, Eq. 
(59) reduces to 

^11 = ^1 + I ^n, Re(a,) + X ^ny M^,) (60) 



= ^1 +^1112^6(^2) + ^112 Mq2) 

+^iii3Re(ai3) + //iii3hn(q3) 

+ ^1122 Re(a22) + //h22 M^2 ) 

+^1123 Re(a23 ) + //n23 Ini(% )■ 



(61) 



If we now consider the nonlinear large-signal ^ 
parameter representation for Z?n, once again assuming a 
two-port network and taking into account the first three 
harmonics, we have 



I^u 



(62) 



/=2,3 



or 



1^1 +^112^2 + = 



% 



+^1211^21 + ^212 ^2 + ^213 ^3 • 



(63) 



Here, ^y^ are functions of all of the harmonics, not just 
the fundamental terms. So for any change in any ap a 
new set of ^y^i will need to be determined. Separating 
the real and imaginary components of the a's, we can 
express eq. (63) as 

/>„ =^„, Re(q, ) + ^„, [Re(^, ) + yTm(^ )] 

+^,„3[Re(a,3) + yTm(q3)] + ^,„ [Re(q, ) + yTm(g, )] 
+^,,,JRe(a,J + yTm(^, )] + ^,,3 [Re(q3 ) + yTm(^ )]. 

(64) 

Once again, the imaginary component of an is omitted, 
with the phase reference such that the phase of an is set 
to zero. 

We can now equate the nonlinear large-signal ^ 
parameters of Eq. (64) to the nonlinear scattering func- 
tions of Eq. (61), with the understanding that this is 
only generally valid for the special case when the non- 
linear large-signal ^parameters are constant for a 
given bias and fundamental drive level, like F^, Gj^y, 
and Hj^y are defined. Normally, however, the nonlinear 
large-signal ^parameters depend upon the higher har- 
monics as well as on the bias and fundamental drive 
level. The implication of this special case will be dis- 
cussed shortly, after Eqs. (61) and (64) are equated. 
Equating the corresponding real and imaginary compo- 
nents of the a wave variables in Eqs. (61) and (64) 
gives 



^11 =^1111 Re(aii) + ^211^1 ■ 



(65) 



/=1,2 

y=2,3 



i=\,2 
y=2,3 
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Additionally, 



!ll2 ~ ^112 ' 7^112 ~ ^1112 ' 
!ll3 ~ ^1113 ' 7^113 ~ ^1113 ' 
!212 ~ ^122 ' J ^212 ~ ^1122 ' 



and 



^1213 ^1123 ' J ^213 ^1123 ' 

Equations (66)-(69) imply 



(66) 
(67) 
(68) 

(69) 



Gun—JHun'Am—Mm ''Qxii —M22 'Q23 —J^n ' 

(70) 

which means 
Re(G„) = Im(//^,) ; Re(//^,) = -Im( G^,). (71) 

Equation (71) satisfies the conditions of the Cauchy- 
Riemann equations [15], 

dlRcjb^)] _ d[lmib^)] 8[Re(^^)] _ ^lm(^)] 
a[Re(a,)] a[Im(a,)] ' 8[Im(c^,)] 5Re(q,)]' 

(72) 

which implies bf,^ must be an analytic function of a^. A 
complex-valued function is said to be analytic on an 
open set W if it has a derivative at every point of W. 
This is generally true only when b/g, is a linear function 
of ay. Thus, equating the nonlinear large-signal ^ 
parameters with the nonlinear scattering functions is 
generally valid only in the small-signal, linear case. 

As we mentioned earlier, Eqs. (65)-(70) are only 
generally valid in the special case when the nonlinear 
large-signal ^parameters are constant for a given bias 
and fundamental drive level, like F^, G^^.^, and Hf^^y are 
defined. Since this is not generally true, the formula- 
tions for nonlinear large-signal ^-parameters and non- 
linear scattering functions are not equivalent. 

We can draw a few important conclusions, however, 
after attempting to equate the two formulations. First, if 
Gig,ij and Hj^^y are allowed to be functions of higher har- 
monics, then only one of them, either Gf^^y or //^,y, or 
equivalently ^y^, is required since Eq. (70) shows that 
they are not independent. Second, if the nonlinear 
large-signal ^parameters are complex constants for a 
given bias and fundamental drive level and are not 



functions of the higher harmonics, the parameters have 
the limitation that they cannot map circles into ellipses, 
but rather can only map circles into circles, as shown in 
Figure 22. This is because S>iji^ is a single, complex con- 
stant rather than a pair of independent complex con- 
stants such as G}^ij and Hi^y. Thus, if ^g^i is not depend- 
ent upon higher harmonics, it acts like a linear ^S-param- 
eter. 
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Fig. 22. If^ijM is a complex constant for a given bias and fundamen- 
tal drive level, it has the limitation that it can only map circles into 
circles. 
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We have shown above that the two formulations are 
not equivalent. Nonlinear large-signal ^-parameters 
are more general than the nonlinear scattering func- 
tions, which are useful in approximating a specific 
class of nonlinearity in a more compact form. 
Nonlinear large-signal ^-parameters have the advan- 
tage of being able to map circles into any arbitrary 
shape, rather than being limited to ellipses. 
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